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L" We discuss the status of the pion distribution amplitude (DA) in con- 

nection with QCD sum rules and experimental data on the 7*7 — > tt° 

£f-\ . transition form factor. Contents: (a) Pion DA in generalized QCD Sum 

Rules (SRs); (b) Light Cone Sum Rules (LCSR) analysis of the CLEO data 

"^J , for the 7*7 — > n° transition form factor; (c) Recent lattice QCD data for 

the pion DA; (d) BaBar data — a challenge for QCD? 

Pion distribution amplitude from QCD sum rules 

The twist-two pion distribution amplitude (DA) parameterizes the matrix 
element of the nonlocal axial current on the light cone [1] 

1 
(0\d(z)^j 5 [z,0]u(0)\n(P)) =iUpJdxe^ z - p ^^- 2 {x^ 2 ). (1) 

Z 2 =0 J 



The gauge- invariance of this DA is ensured by the lightlike gauge link [z, 0], 
inserted between the two separated quark fields. The physical meaning of 



* Presented by the the first author at the International Meeting "Excited QCD" , Jan- 
uary 31-February 5, 2010, Tatranska Lomnica (Slovakia) 

(1) 



A. P. Bakulev et aJ. 



Excited QCD, 2010 



this DA is quite evident: it is the amplitude for the transition tt(P) — > 
u(Px) + d(P(l — x)). It is convenient to represent the pion DA using an 

expansion in terms of the Gegenbauer polynomials C n (2x — 1), which are 
the one-loop eigenfunctions of the ERBL kernel [2J E] ■ This representation 
means that the whole scale dependence in ip 7T (x;fi 2 ) is transformed into the 
scale dependence of the Gegenbauer-coemcients 02 (/i 2 ), a^/j, 2 ), .... 

In order to construct reliable QCD SRs for the pion DA moments, one 
has to take into account the nonlocality of the QCD vacuum condensates — 
as it has been shown in [U [5]. As a concrete example for the nonlo- 
cal condensate (NLC) model, we use here the minimal Gaussian model 

(<j{Q)q( z )) = (qq) e with a single scale parameter A 2 = (A; 2 ), charac- 

terizing the average momentum of quarks in the QCD vacuum. The value 
of A 2 has been estimated in the QCD SR approach and also on the lattice 
00 El E]: A 2 = 0.35 -0.55 GeV 2 . 

The NLC SRs for the (twist-2) pion DA produce a bunch of self-consistent 
two-parameter models at the normalization scale ^ ~ 1.35 GeV 2 : 

^ LC (x; ^) = <p As (x) [l + a 2 (/, 2 ) C 2 3/2 (2x - 1) + a 4 (/i§) C7 3/2 (2* - 1)1 . (2) 



The central point corresponds to 



,BMS 



+0.188, a 



BMS 



-0.130 for 



A 2 



0.4 GeV 2 , whereas other allowed values of the parameters a-i and 
04 — in correspondence with associated values of A 2 — are shown in the left 
panel of Fig. Q] in the form of slanted rectangles [10] . All these solutions 
yield, in accord with (J2J, the same value of the inverse moment of the pion 



DA, namely, 



1 \ bunch 



L\b 



3.17±0.20. This range is in good agreement with 



the estimates derived from a dedicated SR for this moment which can be 
obtained through the basic SR by integrating over x with the corresponding 
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Fig. 1. (a): The allowed values of the parameters a-i and 04 of the pion DA bunch 
@, evaluated at fi 2 = 1.35 GeV 2 for three values of the nonlocality parameter 
A 2 = 0.4,0.5, and 0.6 GeV 2 . (b): Shapes of three characteristic pion DAs — BMS 
(solid line), CZ (dashed line), and the asymptotic DA (dotted line). 
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weight x~ l (at ^ ~ 1.35 GeV 2 ): (aT 1 )^ = 3.30 ±0.30. It is worth empha- 
sizing at this point that the moment (x _1 )^ R can be safely determined only 
with the use of NLC SRs because of the absence of endpoint singularities. 

Comparing the obtained pion DA with the Chernyak-Zhitnitsky (CZ) 
one [TT|, reveals that, although both DAs are two-humped, they are quite 
different, with the BMS DA being strongly endpoint suppressed (see Fig. 
IHa)). 

Qualitatively similar results have been also obtained with the improved 
Gaussian model of the nonlocal QCD vacuum [12j. In that case, the allowed 
region for the parameters 02 and 04 in Fig. HJa) is shifted along the diagonal 
farther to the right with the central point being located near the right corner 
of the allowed (green) region. We emphasize that the BMS model |10| is 
inside the allowed region obtained with the improved QCD vacuum model. 
This means that the characteristic features of the BMS bunch persist for 
the improved bunch as well — in particular the endpoint suppression. 

Analysis of the CLEO data on F 11 * 7T {Q 2 ) and the pion DA 

Many studies [2H [22J 123 121 El El E5] have been performed to determine 
the pion DA characteristics using the high-precision CLEO data |14] on 
the pion-photon transition form factor F 7T - 1 * 1 {Q 2 ). For instance, we have 
used in [15] LCSRs [261 EI] with next-to-leading-order accuracy of QCD 
perturbation theory with the aim to analyze the theoretical uncertainties 
involved in the CLEO-data analysis and extract more reliable estimates for 
the first two coefficients 02 and 04, which parameterize the deviation from 
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Fig. 2. (a): LCSR predictions for Q 2 F 1 * J ^, n (Q 2 ) for the CZ DA (upper dashed 
line), BMS- "bunch" (shaded strip), and the asymptotic DA (lower dashed line) 
in comparison with the CELLO (diamonds, [13]) and the CLEO (triangles, 
[14]) experimental data, evaluated with the twist-4 parameter value S^ w _ 4 = 
0.19 GeV 2 [H] and at /i| Y = 5.76 GeV 2 . (b): The scaled pion form factor calcu- 
lated in the NLC QCD SRs (solid blue line) including nonperturbative uncertainties 
(shaded strip) [IB] . The dark gray strip shows the lattice data [T7| and the exper- 
imental data are taken from |T5| (diamonds) and [T5], [3D] (triangles). 
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the asymptotic expression ip^ s . The upshot of our analysis [151 [27] 1S that 
the CZ DA is excluded at the 4<r-level, whereas the asymptotic DA is off at 
the 3<7-level, while the BMS DA (and most of the BMS bunch) is inside the 
la-error ellipse, and the instanton-based model of Ref. [28] is close to the 
2<r-boundary. Moreover, we found that the CLEO data conform with the 
value of the QCD nonlocality parameter A 2 , = 0.4 GeV 2 . 

Pion form factor and JLab data 

In this context it is worth mentioning our analysis of the pion's electromag- 
netic form factor which employs NLC QCD SRs, Analytic QCD Perturba- 
tion Theory, and NLC-derived pion DAs [16J. The obtained results are in 
excellent agreement with the recent JLab data and also with the lattice data 
of Ref. [T7], as one sees from Fig. E^b), where the light gray strip includes 
the NLC QCD SRs uncertainties, whereas the dark gray strip represents 
the lattice results. 



New lattice data and the pion DA 

Rather recently, new high-precision lattice measurements of the second 
moment of the pion's DA (S, 2 )-k = L (2x — l) 2 ip 7T (x) dx appeared [301 EI]- 
Both groups extracted from their respective simulations values of a% at the 
Schmedding-Yakovlev scale //| Y ~ 0.24, but with different error bars. Re- 
markably, these lattice results are in striking agreement with the estimates 
for 0,2 both from NLC QCD SRs [TO] and also from the CLEO-data analyses 
based on LCSRs [21J.I15J. The improved bunch [12] appears to have an even 
better agreement with the recent lattice results of [31] . 

Confronting NNLO LCSR results for F 7 * 7 -^. 7r with the BaBar data 

In a recent paper by two of us [32], the NNLO QCD radiative corrections, 
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Fig. 3. NNLO/3,, LCSR predictions for Q 2 F 7 - 7 ^ 7r (Q 2 ) for the CZ DA (upper red 
line), BMS- "bunch" (shaded green strip), and the asymptotic DA (lower solid line) 
in comparison with the CELLO, CLEO (the same as in Fig.QJa)) and BaBar [25] 
(brown diamonds with much smaller error-bars) experimental data. Left panel: 
Q 2 < 9 GeV 2 region. Right panel: whole BaBar region of Q 2 . 
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proportional to the /3o-coefncient, have been calculated and included into the 
LCSR analysis of the pion-photon transition form factor. The overall effect 
of these corrections appears to be negative, hence reducing the 7*7 — > it 
form factor by —7% at low Q 2 ~ 2 GeV 2 and by —2.5% at intermediate 
Q 2 > 6 GeV 2 (left panel in Fig. gj). In spite of this reduction, the BMS 
bunch describes rather well all data for Q 2 E [1.5, 9] GeV 2 — including those 
of BaBar |29j . This means that the CLEO (and the low-energy BaBar) 
data are incompatible with wide pion DAs and demand that the endpoints 
x = 0, 1 are stronger suppressed than in the asymptotic DA. 

Surprisingly, the BaBar data |29| contradict this behavior in the high- 
energy region: there they show a significant growth with Q 2 for values above 
~ 10 GeV 2 . This behavior of the BaBar data is clearly in conflict with the 
collinear factorization in perturbative QCD. This is true for any pion DA 
that vanishes at the endpoints x = 0, 1 (see Fig. El — right panel) and the 
explanations in |32j. 

Despite the appearance of several proposals to explain the anomalous 
behavior of the high-Q 2 BaBar data [33J IM1 ESI ES] , we want to emphasize 
that from the point of view of QCD it is not possible to describe them either 
with the inclusion of higher radiative corrections or with higher-twists — both 
types of contributions are negative [151 E21 EI] • 

Conclusions: 

To conclude, the QCD SRs produce an endpoint-suppressed bunch of pion 
DAs in agreement with the CELLO and the CLEO data on i ? 7 * 7 _>. 7r . The 
values of the second moment of these DAs are within the range found in 
recent lattice simulations. These results are also in line with the new BaBar 
data up to ~ 10 GeV 2 but do not reproduce the observed growth above this 
scale. It was shown in [38] that the CELLO and the CLEO data cannot be 
fitted simultaneously with all the BaBar data. 
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